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Abstract: Three models from statistical physics can be analyzed by employing 
space-time determinantal processes: (1) crystal facets, in particular the statistical 
properties of the facet edge, and equivalently tilings of the plane, (2) one-dimensional 
growth processes in the Kardar-Parisi-Zhang universality class and directed last pas- 
sage percolation, (3) random matrices, multi-matrix models, and Dyson's Brownian 
motion. We explain the method and survey results of physical interest. 
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1 Introduction 



Exactly solvable models from Equilibrium Statistical Mechanics in two space-dimen- 
sions are a continuing source of fascinating research. Historically the Onsager solu- 
tion of the two-dimensional Ising model stands out. But the field has progressed. 
The most recent advance is conformal field theory and its probabilistic underpinning 
through the link to the Schramm-Loewner evolution (SLE) [H]. 

Going back to the Ph.D. thesis of Ising, one conventional starting point for an 
exact solution in equilibrium statistical mechanics is to write the Boltzmann weight 
as a product of matrices. Each factor is referred to as transfer matrix. In case of a 
spin system on the two-dimensional lattice Z 2 the matrix elements of the transfer 
matrix are labeled by spin configurations in adjacent columns. Under favorable 
circumstances the transfer matrix can be rewritten as the configuration space kernel 
of exp[— H], where H is bilinear in fermionic creation/annihilation operators labeled 
by the sites of a single column. In such case one says that the model can be solved 
through a mapping to free fermions. The problem of finding the largest eigenvalue of 
the 2 N x 2 N transfer matrix is reduced to the eigenvalue problem of an N x N matrix. 
For the 2D Ising model the free fermion method was discovered by Lieb, Mattis and 
Schultz |4T) . Under less favorable circumstances one can still extract information 
from the transfer matrix by more sophisticated methods, like the Bethe ansatz, the 
Yang-Baxter equations, and the technique of commuting transfer matrices. We refer 
to \62\ [H] for details. For the purpose of our lectures the free fermion method with 
suitable extensions will do. 

To be more concrete, and to anticipate some features to be developed in much 
greater depth further on, let us explain the free fermion method in the context 
of the ANNNI model [46J, to say the 2D anisotropic next nearest neighbor Ising 
model. The spins, with values ±1, are located at the sites of the square lattice Z 2 . 
Because of anisotropy the properties in the 1- and 2-direction are very different. 
It is then convenient to consider the 1-axis as fictitious "time" and the 2-axis as 
"space". Along the time direction the spins interact via the ferromagnetic nearest 
neighbor coupling J , Jo > 0) while in the space direction there is a nearest neighbor 
ferromagnetic coupling Ji and a next nearest neighbor antiferromagnetic coupling 
J 2 - We set J 2 = —Ji/2 < 0. The ground state is then highly degenerate. It consists 
of alternating strips of + spins and — spins parallel to the time direction with a 
width of 2 or larger. 

As the temperature is increased, the domain walls, i.e., the lines separating the 
+ and — domains, become thermally rough. In approximation we postulate that 
a domain wall has up-steps and down-steps of size 1 only and has no overhangs, 
see Figure d^a). To have easier comparison with the models to be discussed below, 
the distance between neighboring domain walls is diminished by one lattice spacing. 
Then we obtain the following line ensemble: There are M domain walls t i— > Xj(t), 
j = 1, . . . , M, with Xj G [-N, . . . , N], t G [-N, . . . , N]. Each domain wall has steps 
at most of size 1, Xj(t + 1) — Xj(t) = 0, ±1 for all j, t, and the domain walls satisfy 
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Figure 1: (a) Domain walls of the ANNNI model and (b) line ensemble of the PNG 
growth model. Note that for the line ensemble (b) there is a top line. Mostly, our 
focus will be on the statistical properties of the top line. 



the non-crossing (non-intersecting) constraint Xj(t) < Xj+x(t). To an admissible 
configuration of domain walls one assigns a Boltzmann weight. A single step of the 
domain wall has weight e~ Jo ^ kBT . Thus the Boltzmann weight for a configuration of 
admissible domain walls is 

exp [ — -j^Tp (number of up- steps and down-steps)] . 

So far we have imposed a hard wall boundary condition in the space and free 
boundary conditions in the time direction. One could also require periodic boundary 
conditions in the t-direction. A further popular choice are in addition chiral bound- 
ary conditions along the space direction. Then the domain walls have a nonzero 
slope on average. 

The low temperature phase diagram of the ANNNI model, in the approximation 
just explained, can by analyzed using free fermions. Since we will have ample oppor- 
tunity to explain the method, no details are needed now. Let me emphasize that, 
while anisotropic, the ANNNI model as defined is translation invariant. Thus the 
interest is in the free energy per unit volume and in the two-point function which 
depends only on the relative distance of the two points. In contrast, the topic of my 
lectures is concerned with line ensembles which are inhomogeneous both in space 
and time. Figure illustrates the difference. Over the past six years, for some in- 
stants even much further back, it has been recognized that there are three physically 
very distinct systems which can be analyzed through inhomogeneous line ensembles, 
namely 

• crystal facets, in particular the statistical properties of the facet edge, and 
equivalently tilings of the plane, 
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• one-dimensional growth processes in the Kardar-Parisi-Zhang universality class 
and directed last passage percolation, 

• random matrices, multi-matrix models, and Dyson's Brownian motion. 

In the first item the description through a line ensemble is rather immediate, as it is 
in the last item when focusing on Dyson's Brownian motion. For growth processes 
the link turns out to be more hidden. 

Our plan is to explain, in fair detail, for the three items listed how the physical 
model is mapped to a line ensemble with determinantal correlations. I want to 
provide the reader with a feeling for the method, in particular its flexibility and its 
limitations. Thereby one has not yet gained any understanding of the statistical 
properties of the model, but one has reached the "trail head" for a hike towards 
the summit by means of rather formidable asymptotic analysis. The details of the 
hike are well documented in the literature and I will make no attempt to duplicate. 
However, I will discuss of what one learns about crystal facets and growth models 
beyond the specific "exactly solved" model. 

The main technical tool will be determinantal point processes, which are dis- 
cussed per se in the following section. The application to statistical mechanics 
models will occupy the remainder of the survey. 

Acknowledgements. My article is based on the Ph.D. thesis of Michael Prahofer 
|Ho] and on the Ph.D. thesis of Patrik Ferrari [11J. I am very grateful for their shar- 
ing of insights and the ongoing collaboration. In addition, I thank Patrik Ferrari for 
supplying the figures. 

2 Determinantal point processes 

The correlation functions of a determinantal point process are computable from a 
single correlation (two-point) kernel. In this respect, although otherwise very differ- 
ent, determinantal point processes are similar to Gaussian random fields. Therefore 
we first recall briefly 

Gaussian processes. In the discrete setting one starts from the family of mean zero 
Gaussian random variables X\, . . . , X^. Their covariance matrix is 

C(i,j) = (X i X j ). (2.1) 

Clearly C(i,j) = C(j,i) and C > as a matrix. Conversely, every such matrix is 
the covariance of a Gaussian process. Higher moments are computable through the 
pairing rule 

2m m 

<n**>= e n^w*),^*)), (2-2) 

1=1 pairings 7r(fc),7r'(fc) fc=l 

where the sum is over all possible pairings of the indices j\, . . . , j2m- An index is 
allowed to appear several times in the list. 
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Since also space-time processes will be in demand, we add the time index t e 
[0,T]. We have then the family of mean zero Gaussian processes Xi(t), . . . ,X N (t). 
As before their covariance matrix is 

C(i,s;j,t) = {X i (s)X j (t)) (2.3) 

considered as a positive operator on ® L 2 ([0, T]). Except for being continuous, 
the time-index t is on the same footing as the space- index j. In particular, the 
pairing rule ()2.2|) is still valid. 

With this background let us discuss determinantal point processes which is done 
in two steps, static and dynamic. 

Determinantal point processes, spatial part. For simplicity let us first assume space 
to be discrete. At each site j G [1, . . . , N] there is an occupation variable rjj = 0, 1, 
where stands for empty and 1 for occupied, which is the reason for the name "point 
process". We prescribe a hermitian matrix R(i,j) such that < R < 1. Then the 
joint distribution of the 77's is defined through the moments 

m 

(H Vjk ) = det (RUtM)^^ , (2.4) 

k=l 

provided the collection of indices {ji, . . . ,j m } has no double points. As for Gaussians 
there is the single correlation kernel R which fixes the full probability distribution. 
However, note that two different R kernels may give rise to the same //-distribution. 
In the applications below, we will meet a real R matrix which is self-similar to a 
symmetric matrix R, i.e., R(i,j) = g(i)R(i, j)g(j)~ l • Clearly, according to (|2.4I) . 
the moments do not depend on whether they are computed from R or R. 

To understand the connection to free fermions let us introduce the fermion alge- 
bra dj, j = 1, . . . , N, satisfying the anticommutation relations 

{a;, aj} = , {a*, a*} = , {a i; a*} = 5^ (2.5) 

with the anticommutator {A, B} = AB + BA. Let 

TV 

// (2.6) 

be a quadratic fermion operator with f) the NxN one-particle Hamiltonian. H = H* 
is equivalent to P) being hermitian, P) = (}*. We then set rjj = a*aj, the fermionic 
occupation variables, and 

m m 

(H Vjh ) = Z-Hr[e- H Ha* k a jh ] , Z = tr[ e -*] . (2.7) 

fc=l k=l 
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It is an exercise in anticommutators to verify that the moments of the occupation 
variables are determinantal with correlation kernel 

R(i t j) = Z-hr[e- a a^a j ] = ((1 + eV% ■ (2-8) 

Later on a prominent quantity will be the probability for the set B C [1, . . . , N] 
to be empty. To compute it one could use ([2.4)1 : however the fermions offer a 
shortcut. Let xb be the indicator function for the set B and let Nb = J2jeB a j a j- 
Then, with P(-) our generic symbol for probability, 

P(77 - = 0,jeB)= lim Z- l ti\e- H e~ XNB ] 

A^oo 

= lim (det(l + e^))" 1 det(l + e^e'^) 

A^oo 

= det(l - XbRxb) ■ (2.9) 

To define determinantal processes over a continuum rather than discrete space is 
as simple as for Gaussian processes [12] • Let rj(x), < x < £, be the point process. 
Every realization is of the form 

n 
3=1 

Here n is arbitrary and n = means no point in [0, £}. Furthermore we prescribe the 
correlation kernel through a linear operator R on L 2 ([0, £]) with R = R*, < R < 1, 
and tri? < oo. Then, provided {xi, . . . , x m } has no double points, one defines 

m 

(l[v(xk)) = det (R(x e ,x e ,)) 1<e/I<m . (2.11) 
k=i 

Extended determinantal point processes. To extend determinantal processes to space- 
time, it is tempting to follow the Gaussian example by introducing a space index j 
and a time index t € [0,T]. Then, according to the definitions above, rjj(t) would 
be concentrated on a collection of points in [1, . . . , N] x [0, T], which is not a line 
ensemble of the form anticipated in the Introduction. Thus, to properly guess the 
correct structure, we turn to the example of n non-intersecting Brownian paths 
pinned at both ends to the origin, which goes back to Karlin and McGregor [2*7] . 
With such a constraint and boundary condition typical paths resemble vaguely a 
watermelon, compare with Figure El More precisely, for j = l,...,n, t i— > Xj(t), 
< t < T, is a one-dimensional Brownian motion pinned such that Xj(0) = 0, 
Xj(T) = 0. In addition we impose the constraint of non- crossing as 

Xj{t) < x j+1 {t) , 0<t<T, j = l,...,n-l. (2.12) 

The corresponding random field over Kx [0,T] is then 

n 

rj(x,t) = J2K x - x 3(t))- (2-!3) 

3=1 
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Figure 2: Pinned Brownian motion conditioned not to intersect, also referred to as 
watermelon ensemble. 



If T and n are of the same order, then the constraint pushes the top line a distance 
T away from the origin, which is to be compared with the typical y/T fluctuations 
for a single Brownian bridge. Since the repulsion originates in a mere constraint in 
the number of configurations, it is known as entropic repulsion. 

Let us first consider a single line, n — 1. If the Gaussian transition probability, 
from x to y in time t, is denoted by 

p t (x,y) = -i=e^ 2 / 2t , (2.14) 

one obtains 

Po,o(a;i(t) e [x,x + dx]) = p T (0, 0)^(0, x)pr-t(x, 0)dx , (2.15) 

where the subscript 0,0 reminds of the pinning to zero at both ends. pr(0,0) _1 is 
the proper normalization. 

Next let us consider two lines, but for simplicity only the point statistics at fixed 
t. By the reflection principle, for arbitrary starting points Xj(0) = Xj, j = 1,2, 
Xi < x 2 , and unconstrained end points, one has the conditional probability 

P(xi(i) e [yi,yi + dy 1 ],x 2 (t) e [y 2 ,y2 + dy 2 ] \ £i(0) = x 1 ,x 2 (0) = x 2 , 

xi{s) < x 2 (s) for < s < t) (2.16) 
= Z(t,x u x 2 y 1 (pt(x 1 ,y 1 )p t (x 2 ,y 2 ) - p t (x u y 2 )p t (x 2 , yi))0(y 2 - y\)dy x dy 2 

with 9(u) = 1 for u > 0, 6(u) = for u < 0, and Z the normalization. Therefore, 
with pinning at and normalizing partition function Z(S), 

^o,o(xi(t) E [x 1 ,x 1 + dxi},x 2 (t) G [x 2 ,x 2 + dx 2 ]) 

= lim Z(5)~ l [p t (0, x^p^S, x 2 ) - p t (0, x 2 )p t (S, x^] 

x [p T - t (xi,0)p T -t(x 2 ,S) - p T _ t (x 2 ,0)p T ^ t (x 1 ,5)]dx 1 dx 2 
= (R(xi,x 1 )R(x 2 ,x 2 ) - R(xi ) x 2 )R(x 2) xi))dxidx 2 , (2.17) 
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where 

R(x, y) = pt(0, x)pr-t(y, 0) (pr(0, O)" 1 + Zj 1 ^) , 

Z 2 = I pt(0,x)x 2 p T -t(x,0)dx. (2.18) 



Clearly, for fixed t and n = 2, the point process rj(x,t) of ()2.13|) is determinantal. 

To extend to general n, still at fixed time t, it is convenient to first introduce 
the Fermi field over M with creation/annihilation operators satisfying the anticom- 
mutation relations 

{a(x), a(x')} = , {a*(x), a*(x')} = , {a(x),a*(x')} — 6(x — x') , i^El. 

(2.19) 

To have a concise notation one also introduces the harmonic oscillator Hamiltonian 
with frequency 1/t as 

fc^-s?^- (2 ' 20) 

and its second quantization 



H t = \Ja*{x){-^- 2+ ^)a{x)dx. (2.21) 



We now proceed as in (j2.17jl . but for n lines. In the limit 5^0 the left factor 
becomes ^ and the right factor ^T-ti where if)^ is the ground state for H t with 
n fermions. Therefore 

m m 

([[vM) = Z T l (4 l ! t \l[a*(x k )a(x k M n) )F (2.22) 

k=l k=l 

with Z T = (v4_ t |#V and (-1-)^ the 

inner product in fermionic Fock space. It 
is an exercise in anticommutators to confirm that the right hand side of ()2.22|) is 
determinantal with correlation kernel 

Rt(x,x') = Z^(^l t \a*(x)a(x')\4 n) }^- (2.23) 

One can use the normalized eigenf unctions of i)t, i-e., the Hermite functions, to 
express R t as 

n-l 

R t (x,x') = J2^Wx') (2.24) 

j=0 

with t = t(T — t)/T. The correlation kernel is the Hermite kernel of order n. 

With these preparations, we are in a position to attempt our goal, namely to 
extend to the point statistics at several times, e.g., to the joint probability distri- 
bution of T](yi,ti) and ^(2/2^2), < t x < t 2 < T. Of course, following the examples 
above, one can work out concrete formulas. They tend to be lengthy and it is more 
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transparent to emphasize the general principle. The n-dependence turns out to be 
simple and we keep n arbitrary, but the reader is invited to verify for n = 2. From 
the limit 5 — > at t = one obtains the ground state ip^ at time ti, and corre- 

(n) 

spondingly from the limit 5 —>■ at T one obtains V'r-t at time T — fa. For the 
propagation from t\ to t 2 , t — t 2 — ti, one has to use (j2.16|) for general n, which can 
be written as the position space kernel of e~ rG restricted to the n-particle subspace, 
where 

If d 2 

G — —- a(x)—a(x)dx. (2.25) 



2 J v ' dx 2 

More explicitly, in position space, 



\Xi , . . . , X n | 6 



tG \y 1 ,...,y n ) = det (p T (xj - Vj'))^,^ ■ (2.26) 



Note that the free particle Hamiltonian G provides the internal propagation, while 
the harmonic oscillator Hamiltonian H t reflects the pinning. Therefore 

(V{yuh)v(y2,t 2 )) = Z T \^l t Ja%y 2 )a(y 2 )e~ TG a*( yi nyiMP)r (2-27) 

withz T = (^i t2 \4: ) )r- 

It is yet another exercise in anticommutators to verify that the expression (J2.27)) 
is determinantal with the correlation kernel 



mT , , - } Z T \^Ja*{x)e-^ G a{x')\^ ) )^ for t > t , 
1 ' ' ' -Zr 1 <4%,|o( a 0e-"-*' |o o*(a;)lW B) >^ tort<f. ^ 



(I2.27J) generalizes to m pairwise disjoint space-time points (yi, ti), . . . , (y m , t m ) which 
are ordered in time, i.e., < ti < t 2 ■ ■ ■ < t m < T, as 

fit 

(tlviVjitj)) = det (/«'(//,• /.r///-/,')) , (2-29) 

3=1 

The moments are still of determinantal form. In contrast to the static rule, for space- 
time points the time order must be respected. On top the extended correlation kernel 
is not symmetric. 

The general structure can be grasped even more clearly by returning to the 
discrete space setting from (|2.5|) above. In addition to the static Hamiltonian (|2.6j) . 
there is the generator 

N 
i,j=l 

for the time propagation. We define 

aj (t) = e tG a je - tG , a*{t) = e tG a* e - tG . (2.31) 
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Then the dynamic extension of the static correlation kernel R from ()2.8j) is given 
through 



J-%[e-"e-^a*(^')] for < t' < t < T , 
' ' J ' ' ^ ■Z- 1 tr[ C - ff e- TO a /f (Oa;(t)] for < t < f < T , 1 ] 

where Z = tr [e^e -7,6 ] . 

In principle, the moments of the corresponding line ensemble are still defined via 
()2.29j) . For general G, one cannot expect that the so defined moments come from 
a probability measure. In fact, while H can be arbitrary, except for H = H*, the 
conditions on G are rather stringent, as will be explained now, where we distinguish 
whether space, resp. time, is either continuous or discrete. 

(1) continuous time, continuous space. In this case g must be the generator of a 
diffusion process, 

g(t) = -d(x, t) 2 -^ - a{x, t)^- + V{x, t) . (2.33) 

As for the watermelon, the line ensemble is constructed from independent lines with 
a weight determined by the propagator generated through g(t) and subsequently 
imposing the nonintersecting constraint. The watermelon ensemble has constant 
diffusion and formally a strong confining potential at t = and t = T, V(x,t) = 
otherwise. As will be explained, the case d = 1, a = 0, V(x) general corresponds to 
the eigenvalues of hermitian multi-matrix models. 

(2) continuous time, discrete space. g(t) is the generator of a continuous time nearest 
neighbor random walk. If space is Z, then 

(ff(t)V)j = -r + (j,t)ip j+1 - r_(j,t)^_! + (r + (j,t) + r_(j,t) + V(j,t))^ (2.34) 

with rates r± > 0. For the polynuclear growth model we will encounter the simple 
random walk, for which r + (j) = r_(j) = | and V(j) = 0. As before, the line 
ensemble is obtained from independent lines by conditioning on non-crossing. 

For discrete time, there seems to be no complete classification. Trivially, one can 
consider the cases (1) and (2) at discrete times t — m only. In addition one finds 
one-sided exponential, resp. geometric jumps. 

(3) time discrete, space continuous. Let us take the ordered points X\ < ... < x n . 
Then in an up-step the new configuration {x'j, j = 1, . . . , n) has to satisfy Xj < x'j < 
Xj+i, j = l,...,n, formally x n+ \ = oo, and the weight is YYj=i e~ 5 ^ x ^~ Xj \ 5 > 0. 
Our example may look artificial, but does turn up in the analysis of the totally 
asymmetric simple exclusion process. 

(4) time discrete, space discrete. An obvious example are nearest neighbor discrete 
random walks. To have a meaningful non-crossing constraint odd and even sublat- 
tices must be properly adjusted, see Example (i) below. In case space is either Z 
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or Z+, the analogue of one-sided exponential jumps are one-sided geometric jumps 
with weight q n , where n is the jump size, n > 0, and < q < 1. One-sided geometric 
jumps will show up for the Ising corner. In this model time i 6 Z, for t < one 
has only up-steps, and for t > only down-steps, while the geometric parameter q 
depends on t. 

From the list above the guiding principle remains somewhat hidden. There is an 
alternative construction by Johansson j2Sl I2H] which avoids fermions altogether and 
is based directly on a determinantal weight for the line ensemble. In the Addendum 
we outline a purely combinatorial scheme, which was devised by Gessel and Viennot 
as based on ideas of Lindstrom. All the examples from our list can be obtained 
through suitable limits of the Lindstrom- Gessel- Viennot scheme, which therefore 
can be regarded as the most general set-up for determinantal line ensembles. 

Having the machinery of determinantal point processes at our disposal, we turn 
to models of statistical physics. They are crystals in thermal equilibrium (Section 
EJ), growth processes (Sections |U|3J) and the eigenvalue statistics of random matrices 
(Section EJ). Physical predictions are extracted from edge scaling. We also indicate 
briefly how through appropriate boundary conditions for the line ensemble further 
cases of physical interest can be handled. 

Addendum: Nonintersecting paths on directed graphs with- 
out loops 

We explain the Lindstrom-Gessel- Viennot theorem in the form stated by Stembridge 
[13] • One starts with a graph (V, E) consisting of vertices V and directed edges E. 
The graph has no loops. A path P is a sequence of consecutive vertices joined by 
directed edges. V(u, v) denotes the set of all paths starting at u G V and ending at 
v £ V. The paths P and P' intersect, if they have a common vertex. Every edge 
carries a weight w(e) and every vertex a weight w(v). The weight of a path P is 
hence given by 

W (P) = JJ w (e) Y[ w(v) (2.35) 
eePnE vePnv 

and we set 

h(u,v)= w ( p )- ( 2 - 36 ) 

P&V{u,v) 

Let us now consider an r-tuple u = {u±, . . . ,u r } of starting points and an r-tuple 
v = {v\, . . . , v r } of end points. Let Vo(u, v) be the set of all non-intersecting r-tuple 
of paths from u to v. u and v have to be compatible, which means that any r-tuple 
of paths in Vq(u, v) necessarily connects Uj to Vj for j = 1, . . . , r. Then the weight 
of Vo(u, v) is given by 

w(P (u,v)) = det (hiu^vj))^.^. (2.37) 
Let us illustrate the Gessel- Viennot scheme by a few examples. 
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(a) 



(b) 



Figure 3: The directed graph for (a) the Aztec diamond and (b) the 3D Ising corner. 



(i) Simple random walks. Here V = Z 2 restricted to the even sublattice and E are 
all nearest neighbor edges, which are then directed either North-East or South-East. 
Their weight is w{v ) = 1, w(e) = w + for e directed NE, and w(e) = w_ for e directed 
SE. 

(ii) Aztec diamond, domino tiling [221121], see Figure a). Here V = Z 2 and E 
consists of all directed edges as in example (i) plus edges of the form v directed to 
v + (1,0). The horizontal edges have weight w , the NE edges weight w + , and the 
SE edges weight 

(iii) 3D Ising corner, lozenges tiling [15] . see Figure 0(b). Here V — (ZxZ)U((Z + 
|) x Z). The horizontal edges are between nearest neighbors, directed East, and 
have weight 1. The vertical edges are nearest neig hbor for (Z + ±) x Z only. If 
t G Z + | is their 1-coordinate, then for r < they are directed North and for r > 
they are directed South with weig ht gH, < q < 1. 

(iv) Discrete time TASEP [21]. The setup is as in example (iii). Only the North 
and South directed edges are alternating with a r-independent weight q G [0, 1). If 
the vertical lattice spacing is e and the weight is q = 1 — 5e, then in the limit e — > 
one obtains the one-sided exponential jumps from item (3) above. 

(v) Six- vertex model at the free fermion point [12]. In the six- vertex model, see 
e.g. [32], one draws only the South and East pointing arrows. Because of the ice 
rule one then obtains a line ensemble. However lines may touch. To achieve non- 
intersecting lines, SE-edges are added but only for the even sublattice. Touching 
is avoided by choosing the SE short cut. There are then three weights, one for E-, 
SE-, and S-directed edges, respectively. Reconstructing the six- vertex weights one 
notices that they satisfy the free fermion condition. Rotating the space-time lattice 
by 7r/4 one arrives at Figure a), only every second horizontal link is missing. They 
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Figure 4: An atom configuration of the 3D Ising corner. 



can be reintroduced, however, at the expense of splitting up the weight. Therefore 
the Aztec diamond is equivalent to the six- vertex model at its free fermion point. 

3 Equilibrium crystal shape 

As a rule crystals in thermal equilibrium are faceted at low temperatures. In this 
section we will discuss a simplified model, which can be analyzed through the method 
of determinantal line ensembles, see [T51 114j. 

We consider the simple cubic lattice Z 3 . Each site can be occupied by at most 
one atom and the occupation variables are denoted by n x — 0, 1. If the nearest 
neighbor binding energy is —J, J > 0, then the total binding energy of all atoms is 
given by 



At zero temperature only configurations of minimal energy are allowed. If exactly 
N 3 atoms are available, then they form a cube of side-length N. For concreteness 
we assume that the cube occupies the sites [0, . . . , N — l] 3 C Z 3 . If only N 3 — M 
atoms are available, with M < N <ti N 3 , then the binding energy is reduced by 
— 3MJ compared to the perfect cube. However, now there are many configurations 
of minimal energy. They can be obtained by successively removing atoms from 
either one of the eight corners under the constraint to cut exactly three bonds in 
each step. Let us focus our attention at the corner touching the origin, compare 
with Figure 0] The atoms missing at that corner can then be enumerated ed by a 
height function h(i,j), i,j > 0, taking only integer values such that 



h{i,j)>0, h(i,j)>h(i + l,j), h(i,j)>h(i,j + l), lim h(i,j)=0. (3.2) 



The atoms occupy the sites {x \ x\ = i,X2 = j,x^ > h(i,j)}. E.g. the perfect cube 




(3.1) 



\x-y\=l 
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corresponds to h(i,j) = for all The number (= volume) of removed atoms is 

V(h) = Y,Hhj). (3.3) 

i,3>0 

In principle one should introduce such a height function for each corner and 
the volume constraint refers jointly to all corners. For simplicity we ignore such 
inessential complications, thus disregard the other corners, and impose the volume 
constraint on h in the form 

V(h) = M. (3.4) 

Every height configuration satisfying ()3.2)1 and ()3.4)1 has the same energy. Hence at 
zero temperature every configuration has the same weight and the model is purely 
entropic. The only task is to count. 

To deal with the volume constraint it is convenient to switch to the grand canon- 
ical ensemble as 

Z- x e- y ^ T . (3.5) 

Here T is the control parameter for the volume, not to be confused with the tem- 
perature. We are interested in a macroscopic volume, which corresponds to large 
T. Then the average volume, average with respect to ()3.5|) . equals (V)t — T 3 and 
the height is typically of order T. Let us replace h by hr in order to remember that 
the height statistics depends on T. We switch from Z 3 to (Z/T) 3 , i.e., to a lattice 
spacing 1/T instead of 1. In the limit T — > oo fluctuations are suppressed and one 
observes a non-random macroscopic crystal shape, in formula 

lim —h T ([uT], [vT]) = h ma (u,v) (3.6) 

T — >oo ± 

with probability one. Here [•] denotes the integer part, u, v > 0, and h ma is the 
macroscopic crystal shape. In Figure El we display a typical atom configuration with 
volume constraint M = 3 x 10 5 . For the true crystal shape of the low temperature 
Ising model with volume constraint one has to imagine a perfect cube rounded at 
each corner as in Figure 03 

To have a nice looking formula for h ma , we choose the volume constraint as 
2£r(3)T 3 with (j? the Riemann zeta function. Let us define 

f(a, b,c) = — I du dv logfa + be iu + ce iv ) . (3.7) 
47r Jo Jo 

Then the set Sq = {(u, v, w) , u, v > , < w < h ma (u, v)} is parametrically given 
through 

Sq = {2(/(a, b, c) - log a , /(a, b, c) - log b , /(a, b, c) - log c) | a, 6, c> 0} . (3.8) 

As expected the equilibrium shape is symmetric relative to the (1,1,1) axis. Let 
V = {(u,v) , e~ ut2 + e~ v ' 2 < 1}. Then 

fcma(u,v) = onR^\D. (3.9) 
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The equilibrium shape has three facets lying in the respective coordinate planes, 
see Figure |S] T> is the domain where h ma is rounded. Near the facet edge, in the 
direction r = v — u, one has 

hUr, t) = | cosh(r/4)7r- 1 2 1 / 4 r 3 / 2 , (3.10) 

valid for small r, where r denotes the distance away from the edge. The 3/2-exponent 
is known as Pokrovsky-Talapov law [3*4] . 

The expression ()3.7|) has a simple physical meaning, for which we switch to 
the coordinate frame with (1,1,1) as 3- axis. As can be seen from Figure 0] the along 
(1,1,1) projected height profile yields a perfect tiling of the plane with lozenges which 
are oriented either with angle 0, or 27r/3, or 47r/3. Conversely, a tiling by lozenges 
such that asymptotically in each of the three segments there is only a single type 
translates back to an admissible height configuration. We now focus our attention on 
a small neighborhood of a point in T>. For large T the curvature can be ignored and 
under projection the tiling is such that the fraction of each type of lozenges remains 
fixed. Again the grand canonical version is easier to control and we assign to the 
three types of lozenges the Boltzmann weights a, b, c, respectively. f(a, b, c) from 
f!3.7|) is the free energy of such a tiling. Note that f(Xa, Xb, Ac) = f(a, b, c) + log A, 
as it has to be. A tiling of the plane corresponds to a flat surface with a non-random 
slope determined through a, b, c and f(a, b, c) takes the role of the surface tension. 

Switching back to the original coordinate frame by this construction one obtains 
the surface tension a depending on the macroscopic slope Vh of the height function. 
The explicit formula is unwieldy and not so instructive. We now give ourselves some 
macroscopic height profile h defined on (M+) 2 . It must satisfy h > and d\h < 0, 
dih < 0. In the limit T — > oo the macroscopic free energy is additive and hence the 
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prescribed profile h has the total free energy 

T{h) = I dxidx2<j{yh{x\,X2)) . (3-H) 

Minimizing T{h) over admissible height profiles and under the constraint of constant 
volume, 

/ dx 1 dx 2 h(x 1 , x 2 ) — 2(ii(3) , (3.12) 

yields the actual height profile h ma as implicitly defined in (j3.8|) . Here (r is the 
Riemann zeta function. 

Our real interest are the small shape fluctuations on top of the macroscopic 
profile. The three facets are perfectly flat, no fluctuations. For the rounded piece 
one can use the standard Einstein fluctuation argument, which means to expand 
f(/i ma + 5h) — ^(/ima) to second order in 5h, 

J(li ma + 5h) — F{h ma ) = / dx dx'V5h(x) ■ Hess cr(V/i ma (a;i, x 2 ))V5/i(x') , 

Jr2 Jul 

(3.13) 

where Hess a is the 2x2 matrix of second derivatives of a with respect to Vh. 
The inverse of the operator appearing in the quadratic form for Sh(x) defines the 
covariance matrix C(x,x'). The assertion is that, for large T, 

h T ([uT],[vT))-Th ma (u,v), u,veV, (3.14) 

become jointly Gaussian with covariance matrix C. In fact, as proved in jjjj, such 
a property holds provided one integrates (|3.14|) against a smooth test function de- 
pending on u, v. Roughly, C is the covariance of a free massless Gaussian field with a 
strength which is modulated by h ma . Note that the fluctuations are only 0(1), thus 
tiny compared to the same number of independent random variables which would 
amount to a size (Vf) 2 . If in flUU) one integrates over a small square in it, v, then 
this spacially averaged height has fluctuations of size logT. 

We have left out the most intriguing fluctuations close to the facet edge. There 
the crystal steps have much more freedom to fluctuate as compared to the steps in 
the disordered zone, which are squeezed by their neighbors. To be able to analyse 
facet edge fluctuations, we have to set up the line ensemble. 

We return to (|3.2|) and use instead of h the gradient lines hi(t), t 6 Z, £ — 
0, — 1, . . ., see Figure Efa). They are defined through 

t = j-i, ht(t) = h(i,j)+e(i,j), (3.15) 

where 

e(i,j)=-(i+j-\i-j\)/2. (3.16) 

Then 

h e (t) < h e (t + 1) , t < , h t (t) > h t (t + 1) , t > , (3.17) 
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with the asymptotic condition 

lim h e (t) =£. (3.18) 

We extend ha to a piecewise constant function on R such that the jumps are at the 
midpoints, i.e., at some point of Z + |. The gradient lines are then non-intersecting, 
in the sense that 

ht-^t) < hi(t) , t el, (3.19) 

compare with Figure Efb). For the £-th line let t^i < . . . < t^^(e) < be the times 
of up-steps with step sizes s^i, . . . , se, k (e) and let < U jk (e)+i < ■■■ < h,k(£)+n(e) 
be the times of down-steps with step sizes — S£ t k(e)+i, ■ ■ ■ > — se,k(e)+n(e)- The volume 
under the height function h is the sum over the "areas of excitation" for each line. 
Dividing these areas horizontally results in 

k(£)+n(£) 

v w= E E s ^M- ( 3 - 2 °) 

£=— oo j=l 

Therefore the Boltzmann weight of a line configuration is 

II eX P["T E S ^M- ( 3 - 21 ) 

^=—00 j'=l 

To prove that (|3.21j) defines a determinantal process we use the directed graph 
from example (iii) of the addendum to Section El The vertices of the graph are 
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(ZU (Z + |)) x Z. The horizontal bonds are between nearest neighbors and directed 
to the right. The vertical bonds are only on (Z + |) x Z. For positive t they are 
directed downwards, for negative t they are directed upwards. To every horizontal 
bond we assign the weight one. To every vertical bond with 1-coordinate t + |, t 
integer, we assign the weight 

ft = g l ^' l ) q = e-V T . (3.22) 

A path on this directed graph has a weight which is given by the product of weights 
for each step. The line ensemble {hi, i = 0, — 1, . . .} is a collection of non-intersecting 
paths on this graph and their weight agrees with (j3.21|) . 

Following the scheme of Section |2] we introduce the variables i](j,t) in such a 
way that 

. 1 1, if there is a line passing through (t, j) G Z 2 , 

V(j,t) = < . ' (3.23) 

I (J otherwise. 

From the construction of the addendum to Section |21 we know that r)(j, t) has deter- 
minantal moments. The covariance kernel follows from ()2.36|) . Let us first consider 
a single up-step with weight q < 1. According to ()2.36|) a single path starting at j 
and being at i one time unit later has the weight 

(t + (g))..= R J **i-j>0, 

|o forz-j<0. V ' 

The up-step transfer matrix T + (q) for a particle configuration at an integer column 
to the next one is then the second quantization of t + (g), i.e., T + (q) restricted to 
the n-particle space equals S a t + (q) (g) . . . ® t + (q)S a , i.e., the anti-symmetrized n-fold 
product. 

The same transfer matrix can be obtained also from direct summation. Initially 
there are n points. They move upwards under the non-crossing constraint. We want 
to compute the Boltzmann weight (y±, . . . , y n \T + (q)\ x±, . . . , x n ) for initial configu- 
ration (xi, . . . ,x n ) = (x) n and final configuration (yi,...,y n ) = (y)n- If {x) n = (y)n, 
i.e., no step at all, one has the contribution 1 for T + (q). If (y) n differs from (x) n by 
a single step, one has the contribution 

-q^a k a* k+1 (3.25) 

fcgZ 

for T + (q), where the minus sign arises from the chosen order of Fermi operators. 
Similarly a difference of two steps results in the contribution 



q2 \ a fci a fe a fc 2 +i<+i • ( 3 - 26 ) 
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Therefore 

T+(q) = ^^J- a ki ■ ■ ■ a k n a* kn +i •••<+!• (3-27) 

n=0 ' ki,...,k„£Z 

Using properties of Schur polynomials the sum can be carried out resulting in 

T + (q) = exp[J2a*Q + (q) iJ a J }, (3.28) 

where 

(s+(q))h = n j —ou - 1) (3.29) 
% j 

with 9{j) = 1 for j > and 8(j) = for j < 0, which is in agreement with the 
previous argument. Note that Q+(q) is not symmetric because of one-sided steps. 

By the same argument, for down-steps only, T_(q) is the second quantization of 
t_(g) = i+(q)* which implies T_(g) = T + (q)*. 

With this result the Boltzmann weight from t to t + 1 is T + (q t ) for t < — 1 and 
T-(qt) = T + (q t )* for t > 0. In the classification of Section |2] the generator of the 
time propagation is time- dependent. 

To obtain the covariance kernel for the point process r)(j, t) two limit procedures 
are still needed. Firstly we let exactly M + 1 lines run from t = — S to t = S 
and require that at ±S the sites [0, —1, ... , — M] are occupied. The line ensemble 
of interest is recovered in the limits M — > oo and S — > oo. The formula for the 
covariance kernel can be found in llj, Eq. (5.39). 

Independent of this specific formula, the line ensemble has a rather striking 
appearance. For large T, the plane is divided into an ordered and disordered zone 
which is bordered by the two lines 

b+(t/T) = -2Tlog(l - e-l*l/ 2T ) , b-^t/T) = 2Tlog(l - \e~^ 2T ) . (3.30) 

In the ordered zone, with large probability, rj(j, t) = for j > 6+ (t/T) + C(T 1 / 3 ) and 
vU't) = 1 for j < 6~ (t/T) — (^(T 1 / 3 ). The width of the transition region between 
ordered and disordered is C(T 1//3 ), but to find out requires a detailed asymptotic 
analysis. 

Let us now focus our attention on a point ([wT], [vT]) of fixed relative location 
inside the disordered zone, i.e., b^u) < v < b^u). For T — ► oo, close to this 
point, the line statistics becomes stationary in space-time. It encodes the statistics 
of the tiling of the plane with lozenges at a fixed relative fraction depending on the 
reference point (u,v) through Vh ma (u, v). On a mesoscopic scale r)(j,t) is averaged 
over regions of linear size T, but still inside the disordered zone. One then recovers 
the Gaussian shape fluctuations for the rounded piece of h ma , as discussed before, 
see ffTTSjl . 

Clearly, the facet edge corresponds to the top line h (t). h (t) has more space 
to fluctuate. Thus its fluctuation behavior is expected to be very different from the 
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lines deep inside the disordered zone. ho(t) is the microscopic edge between ordered 
and disordered. The properly adjusted scaling with T is thus referred to as "edge 
scaling" , which will be explained in Section |S1 

4 Growth models in one dimension: PNG 

For the Ising corner the appropriate non-intersecting line ensemble can be seen by 
inspection. Still, it is a sort of miracle that the physical Boltzmann weight makes the 
line ensemble determinantal. For growth processes the line ensemble is much more 
hidden and to bring it to light is one part of the discoveries over the recent years. 
Of course, the construction works only for very special growth processes. Also, the 
method is restricted to one dimension. Even then it does not yield information on 
temporal correlations. 

To stress the similarity with the Ising corner we consider in this section the 
polynuclear growth (PNG) model in the droplet geometry. We use xeK for physical 
space and T, T > 0, for the growth time, which should not be confused with the 
time for the line ensemble. The PNG model describes the stochastic evolution 
of the height profile h(x,T), which takes integer values only. A point x where 
h(x + e, T) — h(x — e, T) — 1 , e small, is referred to as up-step, while h(x + e, T) — 
h(x—e, T) = — 1 is a down-step. Larger steps do not occur. The height profile evolves 
by two mechanisms. Firstly, up-steps move to the left with velocity —1 and down- 
steps to the right with velocity +1. Physically the idea is that material can easily 
attach once a step is formed. Steps may collide, upon which they simply coalesce. 
Coalescence should be thought of as a damping, or smoothening, mechanism which, 
as in any other nonequilibrium system, has to be counterbalanced by a suitable 
driving force. For PNG it is given through random nucleation events. They have 
Poisson statistics in space-time. At a nucleation event a nearby pair of an up-step 
and a down-step is created, which then move apart according to the deterministic 
rule. 

In the droplet geometry, one imposes initially h(x, 0) = with a single nucleation 
event at (x, t) = (0, 0). The droplet constraint means that nucleation is allowed only 
on the layers with {h > 1}. Clearly, the height will grow faster in the center than 
at the edges x = ±T. In fact, for nucleation intensity 2, one has 

Urn j=h(yT, T) = 2y/l — y 2 , \y\ < 1 , (4.1) 

1 — >oo ± 

with probability one. Thus the macroscopic growth shape is a droplet, which ex- 
plains the name. 

Since h(x, T) e Z, the height profile is determined by the positions of the up- 
and down-steps and it is sometimes convenient to switch to the step world lines. 
We use the relativistic convention according to which the t-axis points upwards. 
Steps have speed 1 (the speed of light). Nucleation events lie in the forward light 
cone {(x,t) | \x\ < t} of the origin only and are Poisson distributed with intensity 



20 



t 




t = T 



x 



Figure 7: Nucleation events, world lines of steps, and the associated height profile. 



2. Each nucleation event is the apex of a forward light cone, corresponding to 
the world lines of the there created pair of an up- and a down-step. The world 
lines annihilate each other upon collision, where the annihilation events have to be 
determined sequentially starting at t — 0. As a result one obtains an ensemble of 
broken lines. They divide the forward light cone {(x,t) | \x\ < £} of the origin into 
layers of constant height. The lowest layer has height 1, since there is a Poisson 
point at (0,0). Crossing the broken line the height increases to 2, etc.. In Figure 
[7|we display an example with flat initial conditions h(x, 0) = 0. The height profile 
h(x, T) at time T records the height along the horizontal line t = T. 

The construction of the line ensemble is based on the observation that at each 
coalescence of two steps one loses information, since there are many ways how a 
particular height profile could have been achieved. To retain the information we 
set ho(x,T) = h(x,T), the PNG profile, and introduce at time T = the extra 
book-keeping heights he(x, 0) = £, I — —1,-2,.... By definition h (x,T) evolves 
according to the PNG rules. In addition, whenever there is a coalescence event at 
line £, it is instantaneously copied as a nucleation event at the same location to the 
lower lying line £ — 1. Random nucleation events occur only at the top line ho. For 
the book-keeping heights he, £ < —1, the steps move deterministically and coalesce 
according to the PNG rule. 

As before the occupation variables rjT(j,t), \t\ <T, j £ Z, are defined through 



The pattern for rjr{J, t) has an appearance rather similar to the Ising corner. There 
is a disordered zone sharply separated from the ordered zone. Above the droplet, 
j > 2\JT 2 — t 2 , in essence, rj T (j,t) = 0, while below, j < —2\JT 2 — t 2 , one has 
VT(j,t) = 1. At the two borders, r] T (j,±T) = 1 for j < and t)t(J, ±T) = for 
j > 0, see Figure E^b). 



VrU,t) 



1 , if there is a line with h#(t, T) = j 
otherwise. 



(4.2) 
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The joint distribution of the line ensemble {he(t,T) , £ G Z_}, T fixed, is deter- 
mined dynamically through the PNG rules. Surprisingly enough, precisely the same 
distribution can be generated also statically. For this purpose let us consider a fam- 
ily {xe(t) , £ G Z_} of independent, time continuous random walks on Z, i.e., xg(t) 
takes values in Z. We require that 2^(±T) = I. The random walks jump to nearest 
neighbor sites only and do so with rate 1. In other words, the right and left jumps 
occur independently at Poisson times with rate 1. Xf(— T) = £ and the random walk 
is constrained to arrive at i at time t = T. These random walks are conditioned 
not to intersect. The conditioned walks are denoted by xe(t). Then hi(t,T) = xi(t) 
jointly in distribution. The proof is not difficult, but requires some notation. We 
refer to |2E] for the details. From the static construction, it is obvious that r] T (j,t) 
has determinantal moments. In fact the correlation kernel has a structure simpler 
than the one for the Ising corner. Let us first consider the point process {rfr(j, 0), 
j G Z}, along the line t = 0. We introduce the one-particle operator 

(fjrVOi = - Vj-i + ^ • (4-3) 

The eigenvalue equation for fj^, fjT^ = \ip^ x \ has eigenvalues A = |,mGZ, and 
eigenvectors = J J _ m (2T), where J n (z) is the Bessel function of integer order n 
PP. The correlation kernel Bt for {t)t{j, 0), j G Z} is given by 

B T (i,j) = Ji-m(2T)J^ m (2T) , (4.4) 

m<0 

also known as discrete Bessel kernel. The distribution of ?7t(j, 0), j G Z, equals 
the positional distribution for the ground state of an ideal Fermi gas on the one- 
dimensional lattice Z with nearest neighbor hopping and a linear potential of slope 
1/T. The first particle is located typically at j = 2T, while the last hole sits near 
to j = -2T. 

The extension to fermionic time uses the one-particle Hamiltonian 

(flVOj = -^j+i ~ *l>j-i + tyj ( 4 - 5 ) 

which, up to the overall minus sign, is the generator for the time-continuous random 
walk xp(t). Then the space-time correlation kernel is 

U ' , ' y ' J |_( e -*fl(i _ i^e^., fort<tf, 1 J 

with \t\, \t'\ < T. 

Already in their seminal paper Kardar, Parisi, and Zhang recognized that 
growth processes can be reformulated as a directed polymer in a random potential, 
which gives the subject an equilibrium statistical mechanics flavor. This suggests 
that the PNG model, hopefully also the associated line ensemble, must have a 
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Figure 8: Optimal leftmost and rightmost directed polymer over Poisson points. 



transcription to directed polymers. The hint comes from the space-time picture of 
the step world lines. According to convention, the space-time diagram is rotated 
by — tt/4. Then the Poisson points u> = {ujj,j = 1,2, . . .} lie then in the positive 
quadrant (M + ) 2 of the plane. The broken lines are parallel to the coordinate axes 
and have roughly a hyperbolic shape. The directed polymer 7 is, so to speak, dual to 
the broken lines. 7 starts at the origin (0, 0) and ends at (u, v) e (K+) 2 . 7 consists 
of consecutive line segments, which have Poisson points (and (u, v)) as their end 
points. 7 is directed in the sense that each line segment must have positive slope. 
In other words, if 7 is at the Poisson point Uj, then the next Poisson point of 7 must 
be in the forward light cone with apex Uj, see Figure |H1 
To each directed polymer 7 one associates the "energy" 

E(j) = number of Poisson points along 7 . (4-7) 

As can be seen from the geometrical construction, the height is just the energy of 
an optimal path. Thus we define 

e(u,v)= max £(7), (4.8) 

7:(0,Q)— 

where the maximum is over all directed polymers from (0,0) to (u,v). In general, 
there are many maximizing paths. From the geometric construction it is obvious 
that 

h(x,T) = e(T-x,T + x), (4.9) 

up to a scale factor \/2 which is compensated by demanding that the Poisson points 
for the directed polymer have density 1. The dynamical PNG model is replaced 
by finding the energy of an optimal directed path in a random potential. If we 
would associate to 7 the Boltzmann weight e /3E ^ 1 \ then the optimization problem 
corresponds to the zero temperature limit (3 — ■> 00. 

Since h$ = h, through the directed polymer we have reconstructed the top line 
of the line ensemble. To find out we recall the space-time picture of nucleation 
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and coalescence events. Let us call u = u/°) the original Poisson points and uo^ 
the corresponding coalescence points. Of course, will no longer be Poisson 
distributed. We now regard as the second generation nucleation events and 
construct from them h-i by the same rules as we did ho from In turn the 

coalescence points cu^ of are regarded as nucleation points for h-2, etc.. If T 
is fixed, eventually no points remain and from thereon he(x,T) = £. 



5 Growth models in one dimension: TASEP 

A second popular growth model is the TASEP (totally asymmetric simple exclusion 
process). Its height function h(j,T), j e Z, T > 0, takes values in 2Z for j even 
and in TL + 1 for j odd. The height differences are 1 in absolute value, \h(j + 
1,T) — h(j,T)\ = 1. In the growth dynamics local minima of h are increased 
independently by two units after an exponentially distributed waiting time. More 
precisely, if j m is a local minimum of h at time T, then h(j m , T) is updated to 
h(j m , T + t w ) = h(j m , T) + 2 with t w the independent waiting time. If thereby a new 
local minimum is created, one assigns to it a further independent waiting time, etc.. 

The name TASEP comes from interpreting the difference r/ T (j) = (l — (h(j + 
1,T) — h(j, T)))/2 as occupation variables, where rj T (j) = refers to site j empty 
and rfr(j) = 1 to site j occupied by a particle. Translating our updating rule, each 
particle jumps to the right after an independent exponentially distributed waiting 
time provided the right neighbor site is empty. "Exclusion" means that there is 
at most one particle per site and "totally asymmetric simple" refers to nearest 
neighbor jumps exclusively to the right. One could modify the model to its partially 
asymmetric version by allowing also jumps to the left, respecting exclusion. In the 
growth interpretation some material would detach from the surface. 

To construct the line ensemble we consider the particular initial condition 

h(],0) = \j\, (5.1) 

which is the analogue of the droplet for PNG. In the course of time the cone fills up 
and 

lim ^h([uT],T) = h ma (u) , (5.2) 

where 

. . I \u\ for \u\ > 1 , 

h ma {u) = t 5.3 
I ±(u z + 1) for |u| < 1 . 

The key for the construction comes from the directed polymer. Let us consider the 
positive quadrant (Z + ) 2 and attach to each site the random variable w(i,j). 

The w(i,j)'s are independent and have a unit exponential distribution. They are 
linked to the waiting times in the growth steps. As before, we introduce a lattice 
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path 7. It starts at (1, 1), ends at (m,n), and at each step it can either move East 
or North. To such an East-North directed path 7 we associate the energy 



E(rf)= J2 w ^3)- (5-4) 

(m)37 

The energy of an optimal path is then 

G(m,n) = max E(^y) . (5-5) 

G(m, n) is related to the TASEP height through 

¥({G(m, n) < T}) = P({m + n < h(m - n, T)}) . (5.6) 

In the spirit of the PNG model one reinterprets G(m, n) as the height of yet 
another growth process h(J, r) by setting 

h(j, t) = G(t - 1 + j, r - 1 - j) , \j\ < r - 1 . (5.7) 

Hence j G Z, r is the discrete growth time, and h(j, r) G R. The growth process is 
defined through the stochastic iteration 

h(j 7 0) = 0, (5.8) 

{max{h(j - 1, r), h(j + 1, r)} 
+ W ((r + j)/2,(r-j)/2), if (-1^ = 1, 
h(i,T), if = -1, 

for |j| < r + 1, 
^(j,T + l) = Ofor \j\>r + l. 



From (|5.2|) one infers that for large r 



2r 

1 +"(77^ 



fe 0'>T) = TTT7732» \j\ ( 5 - 9 ) 



In particular, the height profile has a macroscopic jump of size r at the boundaries. 

As displayed in Figure 01 the dynamics can be visualized by extending h(j, r) to 
a piecewise constant function with steps on the shifted lattice Z + |. In the random 
deposition step the sequence w(i, r + 1 — i), i = 1, . . . , r, is added at every second 
site from left to right to the current height profile h(j,r — 1). In the deterministic 
growth, up-steps move one lattice unit to the left and down-steps to the right. 
Thereby neighboring steps overlap and the corresponding excess mass is deleted. 

This last rule is the door for the extra book-keeping heights hi(j,r). Initially 
hi(j,0) = 0, £ = 0,-1,.... We set h Q (j,r) = h(j,r). Random deposition takes 
place only at the top height. The sidewards growth is carried out simultaneously 
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+ 



r = 1 r = 2 r=3 

Figure 9: A line ensemble for Laguerre growth with shading 'grey' for nucleated 
mass, 'white' for sidewards growth from previous time step, and 'dashed' for excess 
mass from the line above. 



for all height lines and, at the end of the sidewards growth, the excess mass at line 
£ is copied and added at the same location to line £ — 1. In formulas one sets 

ho(j,r) = h(j,r), 

h t (j,0) = 0, (5.10) 

{ht-i(j,r) - h t (j,T) 
+ mm{h i (j - l,T),h t (j + 1,t)}, if (-l) r+J = 1, 
ht-i{j,r), if (-l) r+i = -1, 

for the line labels £ = 0, —1, . . .. 

As for the PNG droplet it is possible to describe the statistics of the collection 
of points {h#(j, t) \£ G Z_, \j\ < r, he(j,r) > 0} directly without recourse to the 
stochastic dynamics as follows. First we have to define admissible point configura- 
tions. Let {xj,j = —n, . . . , 0} be points on [0, oo) ordered as < x_ n < . . . < x . 
We say that {xj,j = —n, . . . , 0} -< {x'^j = —n, . . . , 0} if x < x' , Xj < x'j < Xj + \ 
for j = —n, . . . , —1. Admissible point configurations of the TASEP line ensemble 
have to satisfy 

ht(±T,T) = 0, (5.11) 
{h t (j,T),ieZ-.} -< {htti + l,r),£eZ-}, if |j|<rand(-l) J ' +r = -l, 
{h i (j,r),£eZ_} y Wj + l,r),f6Z_}, if \j\ < rand (-1)^ = 1. 

As with the growth dynamics, the order -< and >- can be visualized by extending 
hi(j, t) to M by setting hi(x, r) = he(j, r) for j — ~ < x < j + |. Then (I5.11|) means 
that the lines he(x, r) do not intersect when considered as lines in the plane, see 
Figure M 

To a given point configuration, alias line ensemble, one associates a weight, which 
is the product of the weights for each single step, which in turn equals e^' 5 '/ 2 where 
5 is the step size. The total weight is normalized to become a probability, which 
then agrees with the probability from the growth dynamics ()5.10|) at time r. 
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To the line ensemble one associates the point process 

vAy,j) = ^2S(hiU,T)-y), y>0, (5.12) 

where j G Z refers to time and y G M+ to space. According to our construction, at 
y = there are an infinite number of points. The point process r) T refers however 
only to points with a strictly positive y coordinate. t] T is determinantal with a 
correlation kernel R T , which is displayed in Proposition 3.3 of JH]. For j = 0, 
t = 2m + 1 , the correlation kernel simplifies to 

R T (0,y;0,y') = K^(y,y') , (5.13) 

where is the Laguerre kernel of order m. In terms of the standard Laguerre 
polynominals L n of order 0, see [T], it is defined through 

m— 1 

Kt(y,y') = L n {y)L n {y')e- y,2 e- y ' 12 . (5.14) 

n=0 

In analogy, the line ensemble {he(j, r) , t G Z_ , j G Z} is called the Laguerre line 
ensemble. 

Not to our surprise we have found again a disordered and an ordered zone sep- 
arated by the line (|5.9|) on the macroscopic scale. Only in this growth model the 
lower border is trivially the line {y = 0}. 

There are two obvious questions. 

(1) Could one choose for w(i,j) a distribution which is different from the exponen- 
tial without losing the determinantal property? One can, but the only admissible 
modification is w(i,j) to have a geometric distribution, F(w(i,j) = n) = (1 — q)q n , 
< q < 1, n G N, compare with example (iv) of the addendum to Section |21 Note 
that thereby one has constructed a family of growth models, mostly referred to as 
discrete time TASEP, which interpolate between PNG and TASEP. In the limit of 
rare events, q — >• 0, the w(i,j) turn into a Poisson process on with constant 
intensity, while for q — > 1, and proper rescaling, the geometric distribution turns 
into the exponential one. 

(2) Assuming that w(i,j) is exponential, is it required that they all have the same 
mean? In fact not, but the determinantal property requires a product structure. 
More concretely, we assume that w(i,j) are independent exponentials with mean 
(w(i,j)) = (ay) -1 . Then it is required that 

tjjj = <jj + bj > . (5.15) 

Before we discussed the special case Oj = ~, bj — ~. Our construction of the line 
ensemble can be repeated in general, only the weights of the line ensemble have to 
be modified. The up-steps are ordered from right to left and the j-th up-step has 
weig ht e-°*l*l, 6 the step size, while the down-steps are ordered from left to right 
with the j-th down-step having weight e - ^' 5 '. 



27 



6 Random matrices and Dyson's Brownian mo- 
tion 



Random matrices is the most unlikely, and at first sight unexpected, item in our list. 
In retrospect the dynamic exponent z = 3/2 of the KPZ equation in one dimension is 
the "same" as the exponent 1/2 for the edge of the density of states according to the 
Wigner semicircle law. A clear evidence for the link was established by Johansson 
|21j . He proved that, in the droplet geometry, the TASEP height above the origin 
has a scaling function given by the Tracy- Widom distribution, which was obtained 
prior by Tracy and Widom [13] as the scaling function for the location of the largest 
eigenvalue of the Gaussian unitary ensemble (GUE). Historically, it was a big riddle 
why the same scaling function appears in such a diverse context. Our resolution 
is on the mathematical side. Largest eigenvalue and height above the origin result 
from the edge scaling of a determinantal space-time process. 

The GUE of random matrices is a Gaussian probability distribution for N x N 
complex hermitian matrices defined through 



Here A = A* is a N x iV complex hermitian matrix. (j6.1|) is understood as a density 
relative to the flat measure dA on the independent coefficients of A, 



and Zn is the normalizing partition function. In convential random matrix theory 
the factor 1/2N in the exponential is taken to be 1. Our units are such that the 
typical spacing between eigenvalues is of order 1, in accordance with our previous ex- 
amples of point processes. Let Ai, . . . , Ajv be the eigenvalues of A. As a consequence 
of (|6.1j) their joint probability density is given by 



Z~ 1 exp[-trA 2 /2N} . 



(6.1) 



N 



dA = \[ dA %i [ [ , 



(6.2) 



i=l l<i<j<N 



N 



|A^(A)| 2 ne" A ? /2Ar 



(6.3) 



with the Vandermonde determinant 



Ajv(A) = det((A i ) J '- 1 ) 1 < lJ < iV 



n ( A i-^)- 



(6.4) 



l<i<j<N 



We regard Ai, . . . , Ajv as point process by setting 



N 




(6.5) 
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7]n(x) is determinantal with correlation kernel given by the Hermite kernel (|2.24J) 
with t = 2N. 

To make contact with line ensembles one has to advance from the static dis- 
tribution (|6.1|) to dynamics. A natural candidate is the linear Langevin equation, 

|^) = -^W + 5W, (6-6) 

where B(t) is N x N complex hermitian matrix- valued Brownian motion. To say, 
B(0) = 0, (B(t)) = 0, and By(t) are complex-valued Gaussian processes with 
Bij(t)* = Bji(t) and independent increments, 

{BijtfBwit)*) = S(t - s)S w S jr , (6.7) 

i, j, = 1, . . . , N. Clearly, if in ()6.6j) A(0) is Gaussian, then so is A{t). 

The stationary distribution for ()6.6|) is the GUE probability measure (|6.1|) . Thus 
a natural choice is to consider the stationary process for ()6.6|) . The eigenvalues 
Ai(i), . . . , Ajv(£) of A(t) never intersect and form a determinantal line ensemble with 
correlation kernel 

f(x\e- tH »P N e* H »\x') iovt>t', 
R N (x,t:x,t) = < lTT x , rr 6.8 

V ' ' ' 7 \-{x\e- tHN {l-P N )e tH »\x') for t<t'. V ; 

Here Hn is the harmonic oscillator Hamiltonian with frequency 1/2N, 

1 / d 1 o 1 \ . 

and P/v is the Hermite kernel, i.e., Pn is the projection onto the first iV eigenstates 
ofH N . 

As in the previous models, the same line ensemble can be constructed statically. 
One starts with N independent Ornstein-Uhlenbeck processes governed by 

^i(t) = --^Vi(t) + b j (t), (6.10) 



j — 1, . . . , N, with {bj(t) , j = 1, . . . , N} a collection of N independent white noises. 
In the time window t G [— r, r] one conditions on the ?/j(t)'s not to intersect. The 
resulting process is denoted by y - (t). Taking the limit r — > oo one arrives at the 

by construction stationary diffusion process {yj (t) , j = 1, . . . , N , t G K}. It is 
indeed determinantal with correlation kernel ()6.8)) . 

From the perspective of the PNG droplet, stationarity looks unnatural. Closer 
to PNG would be the watermelon ensemble from Section El In terms of random 
matrices one sets 

A(t)=B(t)-^B(T), (6.11) 

i.e., each matrix element is a Brownian bridge, in particular A(0) = = A(T). The 
eigenvalues Ai(t), . . . , Ajv(i) of A(t) are determinantal with correlation kernel ()2.28|) . 
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7 Boundary sources 



From the perspective of growth processes the method developed so far has two 
drawbacks. Firstly, while we have rather concise formulas for spatial correlations 
at fixed growth time T, there is no information on correlations in growth time. 
This limitation is intrinsic, since the line ensemble is constructed separately for each 
T. Secondly, one can allow only for very special initial conditions which result in 
surfaces with a nonvanishing macroscopic curvature. While there is some interest, 
for example the Eden growth starting from a single seed builds up an essentially 
circular shape, in most computer simulations the initial condition is a flat surface, 
which then stays flat on average. 

The restriction to curved profiles can be overcome, at least partially through 
the method of boundary sources, which covers several cases of interest. Boundary 
sources can be introduced for PNG, TASEP, and GUE. To avoid repetition we 
explain only the PNG model, which happens to be the most transparent case. More 
details are provided in the recent survey [IS] . The flat initial height has been resolved 
only recently |4L)j . see also f01E|- It is tricky with extra ideas and therefore slightly 
outside this overview. 

We start with the PNG droplet, as explained in Section HI and add additional 
nucleation events at the two borders of the sample, i.e., at x — ±T. The sources 
are Poisson in time with left rate a_ and right rate a+. Clearly, the sources will 
modify the macroscopic shape. But this is not yet on the agenda. Rather, we 
want to understand how the extra sources modify the line ensemble. Switching to 
the directed polymer, the sources generate additional nucleation events on the line 
{v = 0} with the intensity a + and on the line {u = 0} with the intensity a_. 
In the discrete setting, see Section the exponential random variables would be 
modified such that (w(i,l)) = a+, (w(l,j)) = a_, and (w(i,j)) = 1 otherwise. 
Note that this modification respects the product form, if the vectors a, b are altered 
only in their first entry from ^ to a% — al 1 — |, b\ — aZ 1 — \. Hence (w(l, 1)) = 
(ct + + ct_ — a + ct-)/a + a~. Taking the limit of rare events we conclude that the line 
ensemble for the PNG model with boundary sources is still determinantal, provided 
there is an extra nucleation event at (0, 0) with geometric weight of parameter 

A further, physically natural choice would be to place a single source with inten- 
sity /3 at x — 0. In terms of the directed polymer there are now additional Poisson 
points along the diagonal {u = v}, which should be viewed as a random pinning 
potential. For large (3 the directed polymer stays order 1 close to the diagonal. Any 
deviation would be too costly energy-wise. As (3 is decreased there will be longer 
and longer excursions away from the diagonal until the critical point /3 C , when the 
directed polymer depins. It is conjectured that (3 C = [201111, but there are coun- 
terclaims mostly based on numerical simulations of the TASEP jTH] . Unfortunately 
the source at x = is not covered by our methods, since it does not respect the 
product structure. To have a determinantal line ensemble one can allow for a general 
intensity p(u, v) of nucleation events provided it is of the form p(u, v) = p+(u)p-(v), 
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which can be satisfied for the boundary sources but not for the centered source. 

For the PNG model with boundary sources the line ensemble {hi(x, T), £ = 
0,-1,..., \x\ < T} is constructed according to the rules of Section HJ Since the 
source is in operation only for h , one still has /^(±T, T) — £, £ = —1,-2,.... 
If h (±T,T) = n±, then the corresponding weight is (a + ) n+ (a„) n_ . This looks 
diverging for > 1. However the up-steps and down-steps still carry a dx 

volume element. Since they are ordered, one obtains a factor 1/n! in the partition 
function which makes the total weight finite for all a+, a_ > 0. 

We do not provide the details for computing the correlation kernel, see jT^j for 
the TASEP. There is however one element which we want to point out. In the 
fermion formalism one has a product of transfer matrices, e~ tG , and a few number 
operators, like a*(j)a(j), sandwiched between the right and left vectors 
If a + = cn_ = 0, the boundary conditions are he(±T,T) = £ which translate to 
Q+ = Q, = Q, where Q is the state with sites j < occupied and sites j > 
empty. If a + > 0, then only the right end point of the top line ho is lifted upwards. 
Thus the boundary state becomes 

n + = a*(i; + )h, a*(VO=X>;a*0'), ^ = , (7.1) 

jez 

correspondingly for — , where Q is the state with sites j < occupied and sites j > 
empty. For the PNG model the generator G is the second quantization of nearest 
neigbor hopping, which implies that 

e- tG a*(^p + )e tG = e t{Q + + ^V(^+) . (7.2) 

Hence the boundary creation operator can be moved from the border to the number 
operator a*(j)a(j). 

Let us illustrate this simplification by computing the correlation kernel R a+ , a _ 
at t — 0. From Section 0] we know that for a + = = a_ 

Ro,o(j,f) = B T (j,f) (7.3) 
with Bt the Bessel kernel. In general one has to compute expectations of the form 

m 

Z-\^e- TG a{r)\{a*{j k )a{j k )a*^ + )e- TG \n) r , 

k=l 

Z = (U\e- TG a^j-)a*(^ + )e- TG \n)^ . (7.4) 

This results in a determinantal point process with correlation kernel 

Ra+,a-(j,f) = B T (j,f) 

+ (a + a^-\(l - B T )\^ + )y\(l - 5 T )^-),((1 - B T )ij + ) r . (7.5) 
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The boundary sources modify the correlation kernel through a one-dimensional pro- 
jection operator. Thus computationally the resulting difficulties are increased only 
slightly. 

Even without computation one can guess typical configurations of the line en- 
semble. To compute h Q (±T,T) in terms of the directed polymer, it has to reach 
(2T,0), resp. (0,2T). Therefore h (±T,T) ~ a±T. For the line with label -1, just 
below the top line, we need the extra information on how h-i(x, T) translates to the 
directed polymer. It turns out that for h (x,T) + h-i(x,T) one needs to consider 
two directed polymers, both starting at (0, 0) and ending at (x + T, x — T). They 
are required to visit disjoint Poisson points. Then 

ho{x,T) + h^(x,T) = max (#( 7l ) + E{ l2 )) . (7.6) 

71^72 V / 

fi:(0,0)-y(x+T,x-T) 
7 2 :(0,0)-c(a;+T,a;-T) 

An according formula holds for h Q (x,T) + . . . + he(x,T). If x is near ±T, the 
second directed polymer has almost no Poisson points to visit. Thus h-\{x,T) ~ 
2T(1 — (x/T) 2 ) 1 / 2 and the lines with I < —1 form a disordered zone as before. 
If are small, then ho(x,T) will follow closely h-i(x,T) in such a way as 

to join tangentially the droplet. On the other hand for ct + ,a_ large, ho(x,T) ~ 
((a + — a>-)x+(a + + a-)T)/2. Clearly the most intriguing case occurs when h (x,T) 
is still a line segment but touches tangentially the droplet. At the touching point, 
x = x m , h(x m ,T) is expected to have unusual fluctuations. In the picture of the 
directed polymer, it chooses either one of the two boundaries and the fluctuations 
from the boundary portions are comparable in size to the ones coming from the 
bulk. Such fluctuation properties are studied for PNG in arid for the TASEP 
in 

8 Edge scaling 

For growth processes the physical height corresponds to the top line of the line 
ensemble. Similarly, the facet edge of the Ising corner is encoded by the top gradient 
line, see Figure^] Thus our task is to understand the statistical properties of ho. The 
most basic information is the size of typical fluctuations of ho for large T, which 
defines the scaling exponents, and more precisely the scale invariant probability 
distributions for large T, which defines the scaling functions. Of course, the hope 
is that these quantities do not depend on the details of the line ensemble and thus 
are valid for all line ensembles discussed so far. This is not so unlikely, since the top 
line has a lot of space for fluctuations, which tend to wash out microscopic details. 
As guiding example serves a general step random walk, which on a large scale looks 
like Brownian motion with the variance of the step distribution retained as only 
information on the random walk. Even more ambitiously we expect that, e.g., in 
the case of growth models, the scaling exponents and the scaling functions computed 
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here are valid for all growth models in the KPZ universality class. This is in complete 
analogy to critical phenomena, where models fall into distinct universality classes. 
As a rather common feature, concrete computations can be carried out only for one 
specific member of a class. 

The finite T, resp. finite N, line ensemble is determinantal. If we consider T — ► oo 
and focus our attention on a domain close to the edge, the line statistics there must 
be still determinantal. In other words, we only have to study the limit T —>■ oo of 
the correlation kernel with an appropriate scaling of its arguments. Through the 
determinantal property one deduces the limiting probability distributions from the 
limiting correlation kernel. 

To illustrate how the scheme works let us consider the PNG model in the droplet 
geometry. The starting point is the discrete Bessel kernel ()4.4j) which is the projec- 
tion onto all negative energy states of 1)t from ()4.3|) . For simplicity let us study the 
droplet close to x — 0. Then (h(0, T)) = IT for large T and in the line ensemble we 
consider the window j = IT + yT 13 and t = rT a with y, r of order one and a, (3 to 
be determined. Inserting in (|4.4j) and switching to the variable y one arrives at 

(MO GO = -i>(y + t~ p ) - v>(y - r^) + ^(2T + y T^( y ) . (s.i) 



To have a limit one must set 



T 



P = I (8-2) 



and obtains for T — > oo 

T 2 / 3 f)T^(2/) = (-^ + </)#/)■ (8.3) 
Thus under edge scaling T 2//3 f) T goes over to the Airy operator 

The Airy operator has M as spectrum with the Airy function Ai as generalized 
eigenfunctions 

H Ai M(y - A) = XM(y - A) , (8.5) 

see p. In particular the projection onto the eigenstates with negative energies is 
the Airy kernel 

POO 

K Ai (y, y') = / d\Ai(y + A)Ai(y' + A) . (8.6) 
As established with rigor in jHE], one concludes that 

lim T 1 / 3 B T ([2T + yT 1 / 3 ],[2T + y'T 1 / 3 ])=K Ai (y,y') (8.7) 

T — *oo 

pointwise. 
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To have the extended kernel, see (j4.6|) . one needs the scaling limit of e * 9 with 
t = rT a . Since by the argument above the spatial scale is fixed as T 1 / 3 , one infers 

a = | , exp[-tg] = exp[r(d 2 /dy 2 )] . (8.8) 

While the value for a is correct, the complete asymptotic analysis shows that the 
time propagation is governed by the Airy operator, 

lim T l ' 3 B T ([2T + T l '\y - r 2 )] , T 2 ' 3 t ; [2T + T 1/3 (y' - r' 2 )] , T 2 ' 3 t') 

(y\e~ THAi K Ai e T ' HAi \y'} for r > r' , 

.( y | e ^-flAi(i _ ^1)6^^12/') for r < r' 
= K Ai (y,T;y',T>). (8.9) 

The right hand side of ()8.9|) is the extended correlation kernel of a determinantal 
process, which we denote by £(y, r). £(y, r) for fixed r is concentrated on a discrete 
set of points, whose density vanishes as 

(tor)) = ^2/- 1/2 exp[-4y 3 / 2 /3] (8.10) 
for y — > 00 and increases as 

(tfy, r)) ~ — |2/r /2 - -4t cos(4|y| 3 / 2 /3) (8.11) 
7r 47r|y| 

for 7/ — > —00. As a function of r, ^(y, r) is concentrated on non-intersecting contin- 
uous lines, i.e., 



£(y,T)= J] ( 8 - 12 ) 

j=-oo 

with r 1— > ?/j(r) continuous. Since [H Ai , K Ai ] = 0, £(y, r) and the yj(r)'s are stochas- 
tic processes stationary in r. 

The convergence in (|8.9|) to the extended kernel carries over to the convergence 
of the height h(x,T) of the PNG droplet. One infers that 

lim T-^ 3 (h(rT 2 / 3 ,T) - 2T) = y (r) - r 2 . (8.13) 

Since Airy functions are all over, yo(r) is baptized as Airy process [SHj and denoted by 
A(t). Some of its properties will be discussed in Sectional At the moment we recall 
that r refers to physical space and T to growth time. h(0, T) increases linearly and 
has fluctuations of size T 1 / 3 . In the spacial domain of size rT 2 / 3 the height statistics 
is governed by the Airy process plus a systematic downward bending as — r 2 . Since 
the propagator on the right hand side of (J8.9|) and the static kernel are given through 
H A i, the process Air) is stationary. This is physically quite reasonable. In every 



34 



small region of the droplet one has the same fluctuation statistics, provided the local 
curvature (and possibly linear pieces) are properly subtracted. 

For the Ising corner one also obtains the Airy process for edge fluctuations, as 
anticipated. But no simple short cut as for PNG seems to be available. 

It is instructive to repeat the heuristic PNG argument for stationary Dyson's 
Brownian motion (|6.6|) . The confining potential is V(x) = x 2 /2N, which translates 
to the potential U(x) = (x/2N) 2 /2 on the level of the Hamiltonian H^, see (16. 9j) . 
Its first N levels are filled up which yields the Fermi energy Ep = 1/2. The largest 
eigenvalue of Dyson's Brownian motion is determined by balancing potential and 
Fermi energy. Hence 

U(Xi) = Ep , i.e., Ai = 2N , (8.14) 
which is in agreement with the Wigner semicircle law asserting the asymptotic den- 

1 /2 

sity of states as ir~ - (x/2iV) 2 ) , |ac| < 2N. The determinantal process close to 
the edge is governed by the Hamiltonian ()6.9|) linearized at Ai, i.e., by 

1 d 2 1 . 

Scaling as in (|8.1|) . one concludes that (3 = 1/3. The time direction has correlations 
on the scale N 2 ^ 3 . By stationarity of Dyson's Brownian motion, the edge eigenvalues 
are thus governed by ()8.9|) in the scaling limit N — > oo. 

If instead of the potential \x 2 we choose some other potential V(x), the GUE 
generalizes to 

Z^esxp[-NtT(y(A/N))], (8.16) 

where V is taken as even polynomial with positive leading coefficient. As before, 
the distance between eigenvalues is of order 1. ()6.6|) becomes 

±A(t) = ~V'(A/N)+B(t) (8.17) 



and (|6.9|) is modified to 



H N = ~-^ + U N (x/N). (8.18) 



Un has to be chosen such that V = — log^ g , where ip g is the ground state of if at. Un 
depends only weakly on N. For the construction of the determinantal process one 
fills the first levels of Hn, which results in a Fermi energy Ep = 0(1). The edge, 
x e , is determined through U(x e /N) = Ep. If U'(x e ) ^ 0, then the edge statistics is 
governed by the Airy operator. It may happen that U'(x e ) = 0, but U"(x e ) ^ 0, 
say. Then the edge statistics changes and is governed by the Pearcey process, see 
for a detailed study. 
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Figure 10: Probability densities of the Tracy- Widom distribution functions at /3 = 1 
(GOE), (3 = 2 (GUE), and (3 = 4 (GSE). Note that they are not centered. 



9 Universal fluctuations 

Under edge scaling the top line is governed by the Airy process A(r), which so far 
was defined only rather indirectly as Vq{t) through ()8.12|) . We return to have a 
closer look at its properties. Let us first consider a fixed time, say r = 0. Then 

P(y4(0) <s) = P(£0, 0) has no point in (s, oo)) . (9.1) 

Repeating the computation in ()2.9|) . it follows that 

P(^(0) <s) = det(l - P s K Ai P s ) = F GUE (s) (9.2) 

with i^Ai the Airy kernel and P s the projection onto (s, oo). The determinant refers 
to the Hilbert space L 2 (R). It is well defined, since P s K Ai P s is of trace class for every 
s- Fqve is known as Tracy- Widom distribution. The corresponding probability 
distribution is plotted in Figure El Rather than computing the determinant one 
uses that -Fgue(s) is related to the Painleve II differential equation 

u"{s) = 2u(s) 3 + su{s) . (9.3) 

One picks a special solution, the Hastings-McLeod solution, uniquely characterized 
by u(s) < 0. Then F GU e 

(s) = e~ v ^ with 

POO 

V(s) = - / v(x)dx , v(s) = (u(s) 2 + s)u(s) 2 - u'{sf . (9.4) 

J s 

During these lectures we have repeatedly raised the issue of the statistical prop- 
erties of the layer between ordered and disordered zones. This can now be answered. 
The typical width of the border layer is T 1 / 3 . Relative to the macroscopic location 
at one single point the transverse fluctuations are governed by Fque- 
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The Tracy- Widom distribution was first established in the context of random 
matrix theory |44j . It then appeared in a, at first sight purely combinatorial prob- 
lem, namely Ulam's problem, which poses the question to determine the length of 
the longest increasing subsequence of a random permutation [3|, see also [3] for a 
survey. In fact, this problem is identical to the height h(0,T) of the PNG model. 
To understand the connection let us return to the directed polymer, see end of Sec- 
tion HI with starting point (0,0) and end point (T, T). Let uj = (ui, . . . ,u n ) be a 
realization of Poisson points in [0,T] 2 . We label their 1-coordinates, Xj, in increas- 
ing order. Labeling their 2-coordinates, yj, also in increasing order yields u = 
y<r(i) ),■■■, (x n ,Vcr(n))) and thus a permutation (cr(l), . . . , a(n)) of (l,...,n). 
Clearly, n is distributed as (n!) _1 T 2ri e _T2 and for prescribed n every permutation 
has the same probability. For given permutation a we define It as the length of its 
longest increasing subsequence, e.g., the permutation (6, 2, 5, 1, 4, 8, 7, 3) has £t = 3. 
From the geometry of the directed polymer it follows that £t = e(T, T), see (14. 8j) . 
and hence It = h(0,T). Thus objects from random matrix theory made their ap- 
pearance in growths problems first through the height above the origin in the PNG 
model jHHl and independently for the TASEP [21| . both in continuous and discrete 
time. 

More ambitiously, we step to several space points of the PNG droplet for large 
T, which means several fermionic times, T\ < . . . < T m , for the Airy process and 
consider 

P(^(ti) < Sl ,...,A(r m ) < s m ) 

= F(£(xj, Tj) has no points in (sj, oo) , j = 1, . . . , m) . (9.5) 

Let -f^Ai be the extended Airy kernel of ()8.9|) and consider 

9(y - Si)KAi(T h y\Tj,y')0{y - sj) = A ij (y,y') (9.6) 

as a kernel in L 2 (R) <g> C m . Then, again by repeating the computation leading to 
(EH), 

nA(n) <s x ,..., A(r m ) < s m ) = det(l - A) . (9.7) 

Even for m = 2 the determinant in ()9.7|) cannot be computed in the simple form as in 
(19 .3|) and (J9.4)) , see however [21 WJ\ , and to extract information requires considerable 
effort. 

Physically the most robust information is the two-point function 

C(r) = ((A(r)-A(0)) 2 ). (9.8) 

For small r Brownian motion dominates and 

C(t) = 2\t\, forr^O. (9.9) 

For large r one finds a decay as |t| -2 . In |2J 07] partial differential equations for 
the multi-time distributions of (j9.7J) are derived, which can be thought of as a 
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generalization of (j9.3|) . As one consequence 

C(r) = 2a 2 -2|r|~ 2 , r > 1 , (9.10) 
with a2 = 0.813... the variance of .Fgue- 

10 What have we learned? 

One- dimensional growth models in the KPZ universality class. We have recovered 
the dynamical exponent z = 3/2, which comes hardly as a surprise, since it is 
well established through theoretical arguments and Monte-Carlo simulations. Novel 
is the computation of scaling functions along with the insight that they depend 
on the geometry of the growth process [37]. If there is a non-zero curvature on 
the macroscopic scale, the height fluctuations are governed by the GUE Tracy- 
Widom distribution. On the other hand for a macroscopically flat surface, the 
scaling function depends on how flat the surface is prepared initially. The flat 
surface, no fluctuations at all initially, has a scaling function different from a surface 
where initially the height differences are shortly correlated |3*5] I40j . It may happen 
that a flat piece of the surface joins a curved one. The height fluctuations precisely 
at the junction are governed by yet another scaling function. If the surface is semi- 
infinite, bordered by a hard wall, the scaling function changes [T^j. In this way one 
realizes the GOE and GSE random matrix edge scaling distributions, see Figure ITUl 
and many more. 

Facet edge. In the scaling limit the fluctuations of the facet edge are identical to the 
fluctuations in a growth process with rounded profile. The linear size of the facet 
takes here the role of the growth time T. As argued in [Hj, the Ising model with 
volume constraint and at a temperature below roughening should have the same 
facet edge fluctuations. In fact any model with short range interactions and a non- 
zero facet edge curvature is expected to be in the universality class discussed here. 
There are other surface models which are still determinantal and exhibit facets in 
equilibrium |3*3*l 13 Ij . To establish that their fluctuation properties are determined 
by GUE random matrix theory remains as a task for the future. 
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